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Abstract 
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diagonal elements of the anomalous dimension matrix (the diagonal ones are known from deep 
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1 Introduction. 



Conformal invariance |], |2|, [| has been known for a long time as a powerful tool for the derivation 
of the constraints on the form of the Green functions in massless field theory. However, the real 
world could not be conformal invariant even approximately, since otherwise the mass spectrum of 
the elementary particles might not be discrete. When the interaction is switched on the necessity 
to construct well-defined theory forces us to introduce a UV cut-off, which inevitably breaks the 
symmetry. It is reflected in the fact that the field operators acquire the anomalous dimensions and 
the coupling constant becomes running. In the description of the physical quantities new objects — 
local composite operators — appear which define the dynamical content of the former. The com- 
putation of the spectrum of the anomalous dimensions of these composite operators is one of the 
most important problems of the theory, since these quantities determine the scaling behaviour of 
experimental observables, i.e. the moments of the parton distributions or wave functions depend- 
ing on the particular processes in question. These non-perturbative entries satisfy corresponding 
evolution equations, known as Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) [], ^| one 
for inclusive reactions and Efremov-Radyushkin-Brodsky-Lepage (ER-BL) equation |7], || for the 
exclusive ones 0. In order to solve them one has to determine the eigenfunctions and eigenvalues 
of renormalization group kernels. The main problem arises in the study of exclusive processes [[!]] 
where contrary to the forward scattering one needs to handle the whole tower of local operators 
with total derivatives which mix with each other under renormalization. In the leading order (LO) 
this problem was analyzed in a number of papers |7|, || [H| [11 , [12], |14], |15| . It was recognized 



that the operators which diagonalize the anomalous dimension matrix to the lowest order have a 
conformally covariant transformation behaviour at tree level. 

Nowadays it becomes a common wisdom to perform the analyses in the next-to-leading order 
(NLO) of perturbation theory. The self-consistent approach to this problem requires the knowledge 
of the two-loop exclusive evolution kernels. Because of the complications spelled out above, the 
direct calculation of the anomalous dimension matrices or the evolution kernels beyond the leading 
order turns out to be tremendous task and was resolved up to now in the non-singlet sector 
[]16 , |17| , |T8| , [19 , p0[ 2l| ) |22 |. Since the conformal operators diagonalize corresponding ER-BL 



evolution kernels in the leading order of perturbation theory and the off-diagonal part appears 
only in the NLO approximation, the symmetry breaking starts at the order 0(a) and might be 
expected to be evaluated from the Green functions of some quasi-one-loop objects. Therefore, the 
first step toward the accomplishment of the program of the evaluation of the (9(a 2 )-corrections 
to the evolution kernels, is to specify the particular pattern of the conformal symmetry breaking 
in gauge theories to leading order. To this end true conformal Ward identities are required. 
Due to the UV divergences of the theory there exist anomalies, but only part of them could 
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be related to those which appear in the violation of scale invariance ||23||. Note, that this is 



already true in scalar theories, where there is one-to-one correspondence between dilatation and 



special conformal anomalies induced by the renormalization of the fields J24J] only. In gauge 
theories, a priori, one expects in addition a breaking due to the gauge fixing and ghost terms 
in the Lagrangian, but it should not affect the final result provided we are interested in some 
physically relevant quantities. Thus, the limit of non-running coupling constant corresponds to 
the exact scale invariance; however, it is by no means true that the special conformal symmetry 
also respects this limit because there exists a new source of symmetry breaking which can be 
attributed to the renormalization of the product of composite operators and trace anomaly of the 
improved energy-momentum tensor. However, in the study of gauge invariant quantities these 
effects could be always embedded into the redefined conformal representation of the algebra, so 
that the symmetry will be broken explicitly only by the terms proportional to the /^-function of the 
theory. Thus, conformal covariance is restored in the non-trivial fixed point of the renormalization 
group function. This fact gives us the opportunity to apply the powerful machinery of conformal 
symmetry. It has been exploited by us quite recently to give several QCD predictions coming from 
the use of the conformally covariant operator product expansion |25|, [26|, . 

In the present analysis we will construct the formalism and calculate the two-loop results for 
the singlet ER-BL kernels, restricting ourselves to Abelian gauge theory. The main tools for 
these purposes will be the renormalized Ward identities for the dilatation and special conformal 



transformations p8| , [29| , p0| , pT| , p2| and the commutator constraint coming from the algebra of 
the conformal group. The most important feature of the former is that they are derived without 
any assumptions about the possible form of the symmetry breaking and, thus, turn out to be true 
anomalous Ward identities in gauge theory. 

Consequent presentation will be organized as follows. Next section is mainly introductory. We 
summarize there some basic definitions of the main objects of our analysis: singlet conformal op- 
erators, renormalized Lagrangian of the Abelian theory, as well as other generalities related to the 
transformation of the latter under the action of generators of the conformal algebra. In section 3 
we construct renormalized conformal Ward identities for the Green functions with conformal oper- 
ator insertions. In section 4 they are used to derive the commutator constraint on the off-diagonal 
elements of the anomalous dimension matrix. The most important feature of the conformal consis- 
tency relation is that it allows to find corresponding entries in the n-th order from the calculation 
of the, so-called, special conformal anomaly matrix in the (n — l)-loop approximation. This is 
used by us to derive the two-loop results for the anomalous dimensions in terms of the one-loop 
quantities which are evaluated in section 5. These results are also transformed in the language 
of ER-BL evolution kernels and checked in section 6 for the particular set of two-loop diagrams 
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with fermion bubble insertions in the GQ-kernels. Next, we present the reduction formulae which 
allows to reconstruct the diagonal part of the ER-BL kernels starting from DGLAP analogues. 
The final section is devoted to the summary. Trying to make the paper self-consistent as much 
as possible we add two appendices. In the first one we give the Feynman rules used in the main 
body of the paper. In the second one we list the extended ER-BL evolution kernels responsible for 
the evolution of the generalized non-forward distribution functions which attract lots of attention 
recently. 

2 Preliminaries. 

To start with, we review few basic properties of the theory invariant under conformal algebra 
of space-time. The conformal algebra |], [| |3| is a 15-dimensional Lie algebra, isomorphic to 
0(4, 2). Its generators are the Poincare generators V^, M.^ u together with dilatation T> and special 
conformal transformation /C M . The latter could be defined as product /C M = TIV^TZ" 1 , where 1Z 
produces inversion of space-time points x^ — > x^/x 2 . Apart from commutators which manifest 
the scalar and vector nature of T> and /C M , respectively, there are extra non-trivial commutation 
relations 

[V,K lt ]_=ilC lt , [V^V]_ = iV^ [K ll ,P v ]_ = -2i\g ia ,V + M ia ,]. (1) 

The transformations of the field operators under dilatation and special conformal transformations^ 
are expressed by the formulae 

5% = xp, 5 s <j)(x) = % [(f>{x),V}_ = iD(j)(x), ^ 

5^x^ = 2x fl x v - x 2 g liv , S%<f>(x) = i [4>(x), /C M ]_ = iK^x), 

where the parameters of the infinitesimal transformations have been set equal to unity. Using the 
method of induced representation |2j one can construct from the, so-called, stability subgroup the 
differential representation G for the generators T> and K,^ of the conformal algebra. Here we deal 
with the irreducible standard representation, namely: 

D = i(d (j> + x^d^) , K^ = i hd^x^ - x 2 d^ + 2x^x v d v - 2ix v Y,^ , (3) 

where d^ is a scale dimension of the field 0. The particular choice for them will be discussed 
below. The action of the spin operator 2^ on the fermion and vector fields is defined as usual: 

i 

= ^ [7 M > 1v\_ i/>, = KdwAv - gvpA^). (4) 

2 The 5-sign means the variation of the shape of the fields only 8<p — 4>'(x) — (j>(x). 
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In this paper we intend to construct and study the Ward identities which originates from the 
dilatation and special conformal transformations. The simplest way to derive them is to use the 



path integral formalism [^J. |3(J, [IT], [32| and dimensional regularization in d = 4 — 2e dimensions. 
Their generic form looks like 

(5 G [0 3l ]X) + (P 0l ]5 G liS]X) + ([O^X) = 0, (5) 

where <5 G -symbol means some infinitesimal transformation from the algebra and (...) corresponds 
to the path integral averaging with the weight function exp iS. S = J d d xC stands for the action 
defined in terms of the renormalized Lagrangian of Abelian gauge fields coupled to fermions 

C = ~Z 3 {G, v f - i (d.A,) 2 + Z 2 $i $h/> + Z lf j, e g$4i[> + d,ud,uj, (6) 

where G^ v = d fl A u — d u An is the Abelian field strength tensor. Here we have also added a free 
ghost Lagrangian to make L invariant under BRST transformation. In Eq. (0) X = Y\ k <p{xk) is 
a product of elementary fields (entering in the Lagrangian); its variation is expressed in terms of 
the generators Q of the conformal group acting on the field operators: 

k 

5 G X = £ 0(*i)0(x 2 ) ...i [<f>{xi), Q]_... <f>(x k ). (7) 
i=i 

In Eq. (|5]) we have considered Green functions with a conformal operator insertion which is defined 
as two-dimensional vector 

«-{%)■ 

An explicit form of the conformal covariant tree operators is0 





— * < — — y < — (2) (1) 

where d=d + d and D=d — d — 2%[fgA. We have used the following relation t^ u -+G +tJi G l/ J r = 
G+IgQ. = —G+„gQ_. The renormalized operators are expressed through the non- renormalized 
ones in the following way 

[O j i] = Y,Z jh {e 1 g)O kh (10) 

k=0 

where the matrix Zjk(e, g) is triangular and independent on I due to Lorentz invariance of the 
theory. It is given in the minimal subtraction scheme by 



(QQ% QG^\ i 

[gqzggz) with Z = i + -ZM + .... (11) 



3 Throughout the paper we use the conventions adopted by Itzykson and Zuber |33|1 for Dirac 7-matrices and 
Lorentz tensors. 
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Since we are restricted to the leading twist operators ([|) defined on the light-cone, it is sufficient 
to deal consequently with the collinear conformal subgroup SU (1,1) of the whole group which 
simplify considerably the analysis. For that we project the Lorentz indices with two light-like 
vectors n and n*, which are defined along opposite tangents to the light cone, with n 2 = n* 2 = 
and nn* = 1. 

Beyond the tree level we are not limited to any unique choice of the scale dimensions d^, 
rather we can define them as we wish. It turns out, that the most optimal and convenient set is 
da = 1, d$ = §. Note, that the Ward identities are invariant with respect to change of the scale 
dimensions, since the generating functional is invariant under linear transformation of the field 
operators. For the scale dimensions of the fields we have chosen, the variations of the conformal 
operators are very simple^ 

6 s [O fl ] = -(I + 3)[O fl ], 5 C _ [O jt ] =it.{z d(l)Z^} [Ou^l (12) 

where a jk (l) = 5 jk a(j, I) = 5 jk 2{(j + + 2) - (I + l)(l + 2)}. 

In the consequent discussion and all practical calculations we will always employ the light- 
cone position formalism for the evaluation of the renormalization matrices of different composite 
operators. Therefore, they should be replaced for these purposes by the non-local string operators. 
Thus, we introduce the quark and gluon gauge invariant operators 
(Qqv) 




o ^ ( f !/,'.,. h 2 ) c[h-- 2 n) { ■ )il){Kxn) T i>(Kin) { ' }t'[h- 2 ii). (13) 

L^j^i'^) = G + ^K 2 n) jc^/cm), (14) 

where we have omitted for brevity the path ordered exponential in the first entry. 

To simplify the derivation of the Ward identities it is instructive to use the following variations 
of the Lagrangian 

5 s C = -d^x^ + dJC-d^ — 0-(^ + l)^^<V, (15) 
8„£ = dp [{iXpXv - x 2 g lliy ) C) + 2x u d C - d^ — - (d 4 + 1) 



dC , „. dC 



v d{d v <t>) 0(3^) 

They lead to the non-renormalized scale and special conformal variations of the action 

5 S S = ej d d x {O a (x) + O b {x) - Q^{x)} , (17) 

5 C _S = J d d x2x- {e [O a (x) + O b {x) - tf^(x)] + (d - 2)d ll O B (x)} , (18) 
4 We have used the covariant conformal transformation for the gauge field strength tensor given in Ref. Q . 
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expressed in terms of the following equation-of-motion (EOM) operators 

5 S S S S S 

Q G {x) = fy(x) = i>-^ and fy(x) = —if), (19) 

{Qfah — Q$ + and the operators of class A and B 

^ rBRST rBRST 

Oa{x) = ^-{G, v )\ O b (x) = -^-Od^, <T B (x) = -^-QA v . (20) 

Because of the factor e appearing in Eqs. ( 17|Jl8D , the conformal variation of the action can be 
different from zero only due to the UV divergencies. The appearing operators can be renormalized 
along the line of Refs. However, since these are exactly the operators which enter in the 

Lagrangian density, their renormalization could be gained from the study of the differential vertex 
operator insertion g^S and t-j^S. Namely, a little algebra leads to 

[A«] = \O a ] + [O b ] + f2 G , [A«] = ~ [O b ] . (21) 

HereQ 

[A*(x)]=g^A [A«0n)] = ^ and [A ext (x)]^^[C»^(x)]. (22) 

From the fmiteness of the differential vertex operator insertions the counterterms of [Oa(^)] and 
[C_b(x)] can be fixed: 

/ \ i / I <91nZ 3 \ _ , , If d\nZ 2 d\nZ 2 \ n , s 



p B {x)} = O b {x)+^^AQ^{x) 



2° dg J AK ' 2 V dg * d£ 
d\nZ 2 

in complete agreement with general renormalization theorems [35]. These equalities hold up to 
counterterms which are total derivatives of some other operators TZf. By considering the Green 
functions of the differential operator insertion and from the fmiteness of the special conformal 
Ward identity, we can conclude that they are proportional to the gauge fixing parameter £ and 
thus vanish identically in Landau gauge. The operator O b = [O b ] is finite by itself, since there is 
no other dimension three operator with the same quantum numbers which it could mix with due 
to the renormalization. Now the final step can be easily done with the result 

5 s S = _^ [A ^]- (7 [A«] + (7 G -e)^ G + (7^-e)^, 

5 C S = -^[A^-*[Ai) + (d-2)[A^) + ( lG -e)n G - + (^-e)n^, (24) 



5 For simplifying notations, we use the following conventions: for unintegrated operator insertions and EOM 
operators arisen form the conformal variations of the Lagrangian we keep the dependence on the space-time point 
explicit A(x); for the integrated quantities with weight function 1 we just omit this dependence A = J d d xA(x), 
while for weight 2x_ we use A - = J d d x2x-A(x). 
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where the coefficients in front of the operator insertions are the well-known renormalization group 
functions 

fc(g) = ^ = -eg + (3(g), CT = ^ ln £ (25) 
(cr = —27c) and the anomalous dimensions of the elementary fields 

l%b = —fJ>-r- hi Zo, 7g = — M~r- hi Z 3 . (26) 

The first expansion coefficients of ji = ^7i°' ) + • • • and = + . . . in perturbative series are 

A) = ^T F N f , lf = \c F , iP = 2 -T F N f . (27) 

Aiming consequent generalization of the present considerations to the non-Abelian theory we keep 
dependence on the colour factors explicit. In QED we have to set Cp = Tp = 1. 



3 Conformal Ward identities. 

Combining the results, we have obtained so far in Eq. (|5]), we will come to conformal Ward 
identities, which, however, contain unrenormalized products of operators. The aim of this section 
is to renormalize these products [35|. Since the coefficients appearing in the Ward identities 



depend on the regularization parameter e, this renormalization procedure leads to the conformal 
anomalies. 



3.1 Scale Ward identity. 

The renormalization of the dilatation Ward identity is quite straightforward and does not require 
any additional information as one which comes from the consideration of the differential vertex 
insertions. By using the action principle |571 we can get the results for the products [Oji] [A 1 ] from 
the Green functions with the conformal operator insertion [Ojj\: 

i[O fl ][A*\ = i[O jt A g ] - tjs^Z- 1 +z(lj d d xA^x)j^-- - 2P G ) ^} .J^U (28) 

i[O jl ][At\=i[O jl At\ - EU-^- 1 Pki\- (29) 

When substituted into the Ward identity, the terms proportional to the functional derivative 
acting on the conformal operator will cancel upon similar terms coming from the renormalization 
of [0 3l ]Q G : 
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For quark equation of motion we have 



2^{ZP g Z-'} [OJ 



fc=0 



where the projectors on the quark and gluon sectors are defined by 



Pc 



1 




P 



G 




1 



(31) 



(32) 



It should be noted that the definition of the renormalization constant we have used above does 
not correspond to the minimal subtraction of UV divergencies, rather it comes from the mere 
integration by parts in the functional integral for the Green functions. 

Putting all these equations together, we obtain the renormalized dilatation Ward identity 

{[Oji\8 a X) = £ {iz + (/ + 3)1 + 2( 7G - e)ZP G Z~ 1 + 2( 7 ^ - e)ZP Q Z- 1 } ([O kl ]X) 

k=0 3 

+i^-([O jl A g ]X) + iaip^X) - i{ lip - e){[O fl Q^)X) - z( 7g - e)([O if /2 G ]^),(33) 



where 



Iz 



d ry A_l 

-n—zz 

an 



dg 



(34) 



Note, that only the whole anomalous dimension matrix of the composite operators, 7 = 7^ + 
^■IgPg + 2 7 ^-Pq, is gauge independent. Consequently, the off-diagonal elements of 7^ are gauge 
independent too, while its diagonal entries have a very simple dependence on the gauge parameter. 

Taking into account that the improved energy-momentum tensor 0, [39], f|l| in the massless 
gauge theory admits the anomalous trace [[5], |42| 



e w (aO = ^[A*(aO]+a[A«(s)]-(^ (35) 

and inserting it in the above equality ([33]), making use of the relation {[OjiQ^X) = i{[0 j^N^X}) , 
where is a number of 0-particles operator, we come to the familiar Callan-Symanzik equation 
13] for the Green function, albeit with a conformal operator insertion, where the canonical scale 



dimensions of the field operators are shifted by the anomalous ones [B4 



3.2 Special conformal Ward identity. 

The derivation of the renormalized special conformal Ward identity requires new inputs. We define 
the renormalization constants for the product of the renormalized conformal operator and finite 
operator insertions [A l _] as 

i[O fl ][AL] = i[O fl AL] + i E{^}i*PH-i]. (36) 

fc=0 



S 



For the product with EOM operator we have 

i[Oji]^ = i[o jt n;] - 1 i^iz+hkPu-i]- ( 37 ) 

fc=0 

Simple calculations lead to the results 

Z^ = 2Zb(l)P Q Z-\ Z G = 2Zb(l)P G Z- 1 + Z: s , (38) 

where we have introduced a 6-matrix with elements 

b jk (l) = 6 jk {2(l + k + 3)^ - [1 + (-l) j ' k ](2k + 3)}. (39) 

An additional renormalization constant in Eq. (|3~8"D is extracted from the ^-pole part of the 
operator insertion 

i 



J d d x2x. Mx)ja7x) ~ G+p{x) 6G~(x) = 1 ( 4 °) 



^ 8G+PK U I ) /.: I) 

In the non-local representation, to be used below, the left hand side of this equality is given by 



WM = / - GW*)^} (41) 

in terms of the non-local string operators defined above. 

Finally, the renormalized special conformal Ward identity takes the form 

([O^X) = -it {a(l) + T(l)} jk ([Om-i]X) + i^([0^_\X) + iaiP^X) (42) 
-i(d - 2)([O jl A^]X) - - e)([0^\X) - z( lG - e)([O fl ^]X). 

Since we have expressed everything in terms of the renormalized operator products, we can safely 
put the regularization parameter e = 0. Moreover, because the LHS of this equality is finite 
the same property is fulfilled by the RHS. We have defined above the special conformal anomaly 
matrix 



Til) = limli \a{l) - 2( 7G - e)b(l)P G - 2( 7v , - e)b(l)P Q ] Z' 1 



-(la - e)4ff - ^-Z; - 7€ £* - a(l 



-2 



lG P G + ^Pq] b + 2[ZW, 6]_ + 1 + Z e f, (43) 



where the second equality holds because 7 C is finite. 
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4 Commutator constraint and two-loop singlet anomalous 
dimensions. 



Now we are in a position to derive commutator constraint which induces the off-diagonal elements 
of the anomalous dimension matrix. To this end we will act with operator 

d 



D 



H— 1 + 7 

djj, 



on both sides of the special conformal Ward identity (|42|) . This is completely equivalent to 
the use of the commutation relation coming from the algebra of the collinear conformal group 
[D,JC-]_ = iJC- and the combining use of the results (^) and (f42|) . 
For this purpose we need the following commutators 



d d 
d\i dg 

^ d\i <9£ 



g 



d (3 e 



g- 







dg g J dg 



—a) C- 
dg I d£' 



(44) 
(45) 



Since a = —2jc is gauge independent in QED, the RHS of Eq. (fj5|) vanishes identically. One 
should be careful in treatment of these equalities acting on the differential vertex operator inser- 
tions, since they are only valid when both differentials act on the same objects. Naive application 
could lead to an error. An explicit derivation, which makes use of these commutators, gives the 
following renormalization group equations for the operator insertions [A!_] in QED: 



g^ G ) a 



G T (g§f*\ 



(46) 
(47) 



The following steps are quite straightforward. Now we will act repeatedly with the operator D on 
every term in the special conformal Ward identity without taking into account the action of the 
latter on the field monomial X, since at the end these terms would anyway die out by the use of 
the conformal Ward identity. First 



k=0 



k=0 



7,a + 7 c + 2-6 
g 



2(3^- h G P G + ^P Q \b + (3 



d_ r 

dg 



z*m + z 



eS 



(4£ 

([Oia-i]X), 



where we have taken into account that (physical part of) the special conformal anomaly matrix 
is gauge independent quantity. Next 



D^([Ofl^]X) = -(3 ^U) ([OjAi]*) + (3 ^ 7G ) <Psmx) + (3 (Jgy^ {P^X 
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fc=0 



— fr^ 
dg g 



'dg 



jgPg + ^Pq b\ ([0 M -i]X), 



(49) 



jk 



Da([0 J Ai]X)=P QU) {[O fl Ai]X) + a U^li) (P^X) 



k=0 
d 

d/j, 



Dy G ([0 3l Q G }X) = [ ^ lG \ {\o 3l n G ]x) - lG £ lg%-2$ \ (Pu-i\X), 



d 



%Qf*PQb\ ([Ou-i]X), 



k=0 



jk 



Dy^O^jX) = (^*) ([<V£J*> 



d/j 



(50) 
(51) 
(52) 



The result of Eq. ([51]) manifest the non-minimal character of the subtraction, we have mentioned 
previously. Since the LHS of each equation is finite we have omitted the divergent contributions 
on the RHS. All other terms in the special conformal Ward identity vanish by the action of the 
.D-operator. Owing to linear independence of the Gegenbauer polynomials we can safely omit the 
Green function which enters in the special conformal Ward identity and get the constraint 



^a(0+7°(0 + 2-S(/) 
g 



0. 



(53) 



This equation shows that the off-diagonal part of y in n-loops can be obtained from the knowledge 
of the (n — l)-order of the special conformal anomaly matrix y c and the /3-function. For instance, 
in the leading order we should keep only the a-matrix in Eq. (|53D as compared to the last two 
entries since their perturbative expansion starts at 0(a). Therefore, from the diagonality of the 
a(/)-matrix it follows immediately that all entries in (the first coefficient in the expansion 
7 = -^y^ + (^) y^ + • • •) are diagonal too, and coincide with the moments of the DGLAP 
evolution kernels. 

In the NLO approximation the constraint (|53"D fixes unambiguously the off-diagonal elements 
of the anomalous dimension matrices: 



QQ • ND(1) 
Ijk 

QG ND(1) 
Ijk 

GO ND(1) 
Ijk 

GG t ND(1) 
Ijk 



a(j, k) 
1 

a(j, k) 
1 

a(j, k) 
1 

a(j, k) 



Ijk + PoOjk I + 7 j Ijk Ijk Ik 



QQ^QG^Q) _ QQ^c(0) QG ^(0) 



GQj°) ( QQC(0) 
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gq<°)qq(°) , gg(v)gq<o) 



Ijk 



GQJ0)QGC(0) _ GQC(0) QG (0) 
ij Ijk Ijk Ik 



(54) 
(55) 
(56) 
(57) 



Adding the two-loop results for the anomalous dimensions y^ of the local operators without 
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Figure 1: Generic form of the one- loop diagrams for the renormalization constants Z*. Grey 
blobs correspond to the usual (Eq. ( |A.5| )) or modified (Eq. ( |A.4j ) ) Feynman rules depending on 
the quantity in question. The crossed circle stands for the usual non-local string operator or for 
the modified operator insertion defined by Eq. ([|l|) which determines the w-matrix. 

total derivatives^, which are available for about two decades from the calculations of unpolarized 
[f44f and, recently, from polarized |I5| deep inelastic scattering splitting functions, we obtain the 
complete entry 

which governs the evolution of composite operators in Abelian gauge theory to NLO. 



5 Evaluation of the conformal anomaly matrix. 
5.1 Local anomalous dimensions. 

In this section we evaluate all one-loop renormalization constants which enter in the special con- 
formal anomaly matrix and determine the anomalous dimensions AB 7^ C D< ' 1 ' ) • 

First, in the leading order 2Z^ = f-'Jz Second, the calculation of the diagrams depicted in 
Fig. |I] and the extraction of the pole part from the renormalization constant Z* = \Z^ + ■ • • 
gives: 

™zf = ~* B ifk io,A + B. (59) 



Here we took into account that according to the Feynman rules ( |A.4| ) these diagrams are given 



by the derivatives with respect to the external momenta of the familiar one-loop graphs; and that 



6 Due to the particular normalization of the Gegenbauer polynomials the QG and GQ results of Refs. [ p4[ 
should be multiplied by appropriate factors, i.e. 6/j and j/6, respectively. 
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the differentiation of the Gegenbauer polynomials leads to the 6-matrix: 

G*b + db) + ^ (tit) = £ + M-q (|^) ■ (60, 

Defining also the one- loop approximation of the constant Z* s = -Z*g + • • • in terms of the 
w-matrix, 

AB zi 1] = y, ab ^ (61) 

we finally obtain the result for the special conformal anomaly to leading order: 

Ijk = -bjh Ik + w jk . (62) 

The first entries, which are the usual anomalous dimensions of conformal operators, are well known 
and read 

QQ ^T = -^(3+ (j + 1) 2 (j + 2) ' 4 ^ + 2)+4 ^ (1) ) (63) 

qg (o) -MN/Tp ij 2 + 3j + 4, for even parity 

7j " i0' + l)0' + 2)(j + 3) \j(j+3), for odd parity 1 j 

gq (o) = -<^f [j 2 + 3j + 4, for even parity 

7j " 3(j + 1)(j + 2) \ j(j +3), for odd parity 1 ' 

GG lf = 2 7 g> = (66) 
The calculation of the w-matrix for the QQ-sector can be found in Ref. [p2| . For the GQ-channel 



we employ the non-local representation for the evaluation of the corresponding kernel. Simple 
calculation provides the result 

W r ( Kl ,K 2 ) = ~ C dy \ V dz GQ )C»{y,z) Q O r {yK X +y^zK X + zK 2 ), (67) 
e 27r jo jo 

with the same kernel in the parity even and odd channels 

GQ JC w (y,z) = C F {6(z)-S(y)}- (68) 

Note that this kernel is gauge independent as it should be (the ^-dependent piece we have found 
in the calculation is symmetrical with respect to the permutation y <-> z, and vanishes because 
of the symmetry properties of the operators). The evaluation of the conformal moments of this 
simple kernel is quite straightforward and we get 

(2k + 3)a(j, k) 



(k + l)(k + 2) 
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X 



A jk - + 2) + ^(1) + 4(k + l){k + 2)-%- 

a{j, K) 



with A jk = v ( 3 + k 2 +A ^ - i> {^r) + 2 ^ - fc ) - w + 2 ) - 



(69) 



1 + (-l) J '~ fc 



1 (2k + 3)a(j, jfe) 
6 (A; + !)(* + 2) 



In addition ® G w = GG w = 0. Denning new matrices djp. = bj k /a(j,k) and 
can rewrite the off-diagonal elements: 

«vr> = ( QQ >r - {a - q v} + °%) - - ^r) v*^ 

™ mi™ 



(70) 

ffifc = w jh /a(j,k), we 



QG ND(1) 

Tjfc 

GQND(l) 
GG ND(1) 

7jfe 



= - ( QQ yf - 9a yf ) ^ qg t1 0) - Q % fc QG 7i 0) , 



= ^Q 7 (0) _ GQ^(O) 



K;o)_ G q 7 (o) 



3 7i 0) ) d jk {ft - ^yf} + ^(0)QQ p . fe + GQ, 
7fc J "jfc 7fc yj/k 7fc • 



, 1 '"'/,/,• (71) 
(72) 

^ + G % fc { /3 °- Q ^° ) }' (73) 

(74) 



5.2 Evolution kernels. 



The representations (ffT|-74) allow us to write the off-diagonal part of the evolution kernels as 
convolution of known kernels. The LO kernels AB V^> in the expansion V(x,y\a) = ^-V^°\x,y) + 
(^r) V^~\ x iy) + • • • are diagonal with respect to Gegenbauer polynomials and their eigenvalues 



are - AB 7 f/2: 



1 / C 2 \ 
dx\ i \(2x 

V cu J 



,(°) qgAo) ■ 

3 > 3 



'QQy(O) QGy(0)\ W 
GQy(0) GGy(O) J( X 'y) = - 2 ^GQ^O) GGy(°) 



3 

(7/ 

5 



(2y - 1). (75) 



These results are already known or can be derived from the DGLAP kernels by the procedure 
described below. In the restricted region \x,y\ < 1 these kernels have the general structure 

A = B 



abtM 



r u >(x,y) = 6(y-x) AB F(x,y)± 



y^y 



for 



A^B 



with 



QG F(x,y) 
GQ F(x,y) 
GG F(x,y) 



x 



y 



\ 2x — y — 1, for even parity 
2/^2/ t f° r °dd parity 

x ( 2y — x, for even parity 



2AT/ZV — 



for odd parity 



-y<*(z-2/)- 



(76) 

(77) 
(78) 
(79) 
(80) 
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Here one comment is in order to clarify the differences of the even parity QG and GQ kernels 
available in the literature. Indeed these kernels are not uniquely defined and different expressions 



could be obtained from various calculation techniques [11, 12, 47, 48]. However, apparently 



different functional dependence on the momentum fraction x, y leads to the same eigenvalues in 
the physical sectors of the corresponding channels. The only difference arises in the space of 
the non-physical moments. Namely, the kernels derived by Chase JTTJ] (and listed above) possess 
the off-diagonal matrix elements in the basis of the Gegenbauer polynomials in the non-physical 
sector, while the kernels derived below 

.,2 



^ Fimpr(%-i y) 



2N f T F ^-(-l + 2x-Axy), 



GQ F impi (x,y) = C F — (-l + 2y + Axy) (81) 



V 



are exactly diagonal. The difference can be attributed to the following additional pieces A® F(x, y) 
= 4N f T F ^(2y - 1) and A GQ F{x,y) = 2C F xx(2x - 1) which should be added to the results of 
Chase to get the equations given behind. Note also that the kernels we have listed (taken without 
colour factors) respect the SUSY relation p9| : 

GQ V i0 \x, y)yy = QG V ( °\y, x)(xxf. (82) 

Now we are able to reconstruct the off-diagonal parts of the singlet two-loop ER-BL kernels 
from the local anomalous dimensions we have derived above, i.e. and ( AB "ff^ — AB 7^)djk 
matrices. In the QQ-channel the ^-kernel is given in Ref. P2"| , while for the GQ-channel we have 
been able to perform the summation of the Gegenbauer polynomials by the procedure discussed 
below in section [/]: 



GQ, 



V(x,y) 



-C t 



-c f 



6(y - x) 



x 



X 



6(y-x)-(2x 

y 



J X 

\v 

x - 

y - 



> X 

x I 

y\ 



83) 



34) 



Intuitively, one could anticipate this last result. For that one should note that in Eq. ( [70|) the 
factor l/(k + l)(k + 2) gives just the eigenvalues of the ER-BL kernel in scalar 0jL -theory. Next, 
the off-diagonal matrix elements of the b- matrix are (2k + 3). According to Eq. (|60D, the latter can 
be induced by derivative of the Gegenbauer polynomials C%(2x — 1) with respect to the argument 
x times the factor (2x — 1). Since the derivative increase the index and reduce the order by one 
unity, we come to the equation we have just presented. 

Note that the kernels ( |83"| , |84"D possess also diagonal conformal moments. They can be subtracted 
by applying the projection operator (X—D), where T> extracts the diagonal part of any test function 
r(x,y): 



dxC u A2x 




r(x,y) 



k=0 



1 



Cl(2y-l). 
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The kernel AB g(x,y) that possesses the ^g^-moments is defined as follows 

- ' 2 ~ 

f l dx( C I )(2x-l)(l-V) AB g(x,y)= J j2 AB gJ °! W-l). (85) 

The factor ( AB 7j° — AS 7^ )djk corresponds to the expected breaking of the diagonality of the 
NLO exclusive kernels [H| and it is related to the one- loop renormalization of the composite 
operators. These contributions can be embedded into a redefined conformal representation, which 
differs from the original one by the shift of the scale dimensions of local operators by the anomalous 
ones jj. This means that new eigenfunctions should be orthogonal to each other with the weight 
factor {yy) n+lj i with n — 1,2 corresponding to quarks and gluons, respectively. This gives us a 
guide how to construct new kernels V$(x,y) from the leading order ones. For that we take into 
account that Vs(x,y) should be symmetrical after multiplication by the factor (yy) n+s : 

AA V s (x,y)(yy) n+s = AA V 5 (y, x)(xx) n+s , (86) 

with n = 1, 2 for A = Q, G-channels, respectively, and 

GQ V 5 (x,y)(yy) 1+s = ^ G V s (y, x)(xx) 2+s , (87) 

for mixed ones. These equations tell us that the two-loop off-diagonal terms V(x, y) should be 
proportional to the first term in the Taylor expansion of Vg(y, x) with respect to 5. This recipe 
leads to the following dotted kernels 



ABf, 



V(x,y) = 9{y-x) 



AB F(x,y) ln- + A AB F{x,y) 

y 



(x-*x\ . [A = B 
± U-y ' for 



Note, however, that simple logarithmic modification of the leading order result does not work 
for the parity even mixing channels. The main difference from the other cases comes from the 
fact that Chase's kernels possess off-diagonal elements in the non-physical space of conformal 
moments, which immediately manifest themselves when these kernels are multiplied by factor 
ln(x/y). However, the above symmetry consideration could not fix unambiguously the modification 
of extra pieces AF and, therefore, it requires an explicit analysis. We find the following results 
for A AB F: 

A GQ F V = 2C F [xlny-xlnx], A QG F V = 4T F N f — [xlny-xlnx]. (89) 

Combining the equations (|83| , p4| , |88| , p9D , the off-diagonal parts of the NLO evolution kernels 
read: 

qq v ND(i) = _ (J _ (QQy g, foQyW + + QQ g QQyW _ QQy(°) QQ g 
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Figure 2: Diagrams giving rise to the (3q dependence in the NLO evolution GQ-kernel. 

+ QGy g, GQy(O) _ QGy(°) ^ GQ g ^ ^ (9Q) 

Q Gy ND(l) = (QQy ^QGyiO) + QQ g ^QGy(°)^ (gl) 

GQ T/ (0) 



GQ^ND(l) = _(2_ v ^GQy ^^QQy(°) + f^j^ _GQy( 



(92) 



GGy ND(l) = _ ( J_ P) ^QygQG F (0) + ^gQ Gy (0)^ (g3) 

Here we use a shorthand ® = Jg 1 ote. 



6 Fermion bubble insertions in the singlet kernels. 

In this section we present the check of the self-consistency of our approach by explicit evaluation 
of the fermion vacuum insertions in the GQ-evolution kernels. From the conformal constraints 
we know that corresponding contribution to the off-diagonal part of the evolution kernel is of the 
form 

A G Q y^a) = _|(j _ V ) [ G ^V(x, y) + a %x, „)) . (94) 

Therefore, evaluating the diagrams in Fig. [| we should immediately obtain the desired form of 
the dotted and g- kernels. To this end we have exploited the formalism developed in Refs. [jlS , pOfl , 



which is the most appropriate for the multi-loop calculations of the evolution kernels and makes 
use of the, so-called, a-representation for the Feynman propagators. Corresponding calculations 
run along the same line as for the one-loop diagrams but now we insert the fermion vacuum 
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polarization ili^k) into the gluon lines (see Fig. |2|) 

(*2 \ e 
3^2) (ffg^-knK), with n(e) 



la 6r(l + 6)r 2 (2-e) 



(with /i* 2 = 47T/i 2 ), which results to the mere shift of the power of the denominator of the gluon 
propagator by the magnitude e. Since we are interested only in the off-diagonal part of the 
corresponding diagrams we can safely omit, after subtraction of sub-divergences, the contribution 
coming from the surviving overall constants, i.e. 11(e) etc., since it would be proportional to the 
leading order evolution kernel. The remainder of the graphs leads to the integrals 



e=0 



J dai da 2 [a{ + a e 2 ] j Vi(a u a 2 ) + v 2 (a 1 , "2)7^ | S(x - ya x - ya 2 



(96) 



where the functions Vi read for the parity even and odd channels 



V f 

v 1 [ai, a 2 



vf(a l} a 2 



= -a 2 + y{a 2 + ax), 

= -a 2 + a 2 + 2y(a 2 - ai), 

v 2 y (ax,a 2 ) = v 2 {ai,a 2 ) 

= a 2 a 2 + y (ax + a 2 - 2a 2 («i + oc 2 )) - y 2 [a\ + a\ - 2axa 2 



To perform the final integrations we need the following formulae [50|: 



da5 (x — ya) = Q^i ( x > x ~ V)i 



(97) 
(98) 

(99) 
(100) 



for the first one, and 



f 1 1 f 1 

/ daaa T B 1 ) 1 (x—ya, a—x) = < - 

Jo t + 1 \y 











r +1 - 











x T+] 



9(y — x) H —Q[x — y) 



(101) 



for the second. Putting e = in the square brackets in Eq. (|96|) , we immediately obtain twice the 
leading order kernels^. Evaluating the integral for arbitrary e and differentiating it with respect 
to e, we finally obtain the contribution to the two-loop kernel 

A) 



A GQ V w (x,y) = -fj^F(x,y)+ GQ F(x,y) \nx- GQ F(x,y) In 

+ [ GQ F(x,y) + GQ F(x,y)] In j 1 - 



x 

y, 



(102) 



6(y - x) 



x — > x l 



with 



x \ 1 - 2y, for even parity 
F[x,y) = 2C F - \ 

y I 1 — Ix, tor odd parity 



(103) 



To be precise, in the parity even channel, we get Chase's kernel. 
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while the other functions are given by Eqs. ( |T8| , |?9"D . Obviously, these expressions consist of both 
diagonal and off-diagonal parts of the evolution kernel. Subtracting the sum of the dotted and the 
^-kernels derived in the preceding section it is a trivial task to convince oneself that the remainder 
is diagonal in the basis of the Gegenbauer polynomials. Namely, this provides a new diagonal 
element 



U D (x,y) = { GQ g(x,y) + GQ F A (x, y) lny - GQ F A (x, y) In 



x 



(104) 



+ 



GQ F A (x, y) + GQ F A (x, y)] In M - - j }% - X )-\ X ^ X _ 



and an equivalent representation for the ^-kernel in the GQ-channel 



GQ 



v{x,y) 



-C f 



9(y -x)ln 1 



x 



y , 



x — > X 

y^y 



(105) 



which holds true for the off-diagonal matrix elements. Note that it is most close to the functional 
form of ®Qg. 

Finally, we can conclude that an explicit calculation supports our results derived from the use 
of the conformal consistency relation (|53|). 



7 DGLAP ER-BL reduction. 

Aiming an explicit numerical integration of the two-loop evolution equations one is interested in 
the momentum fraction splitting kernels. While the off-diagonal parts of these kernels are known 
explicitly, it is not the case for the diagonal entries. In this section we outline the formalism which 
fills this gap. We have been able to restore the one-loop ER-BL evolution kernels analytically. The 
reduction formulae we have derived here are sufficient to find numerically the diagonal part (in 
the basis of the Gegenbauer polynomials) of the exclusive singlet kernels from the corresponding 
DGLAP analogues. 



7.1 QQ channel. 

The formal expression for the diagonal part of the ER-BL evolution kernel can be given by the 
following formula (w(x\u) = (xx)"~2) 

1 00 m(x\u) 

V D (x, 2/) = --E -wU C K 2x - WC]i2y - 1), (106) 

where Nj(v) = 2~ iu+1 T r2(u)(v+f)j\ ^ s a norma hzation factor and jj are the anomalous dimensions 
of the conformal operators (see, for instance, Eqs. ( p3| - |66| )). 
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Figure 3: The integration contour for the evaluation of the integral ( |1U9| ). There is a n-th order 
pole in the point u = b and the cut on the real axis from —1 to 1. 

Therefore, the main problem is to perform the infinite summation and construct the generating 
function for the exclusive evolution kernels. Using the Gegenbauer's summation theorem for CJ- 



polynomials [51] this goal can be easily achieved. We have evaluated the following conversion 



function from the forward to non-forward limits: 



2 ((2x - l)(2y - 1) - A^¥yfj z + z 7 
'u + l,u 



-v-l 



\ 2v 



lfty/xxyyz 



2 ((2x - l)(2y - 1) - ^Jxxyy) z + z 2 



In order to find the ER-BL evolution kernel we have to convolute the above equality with the 
DGLAP splitting function and evaluate the integral on the RHS of this equation. For this purpose 
we use the Euler's integral representation of the hypergeometric function, so that the above formula 
looks more compact (here a(u) = (2x — l)(2y — 1) + 4uy/xxyy): 

S ^m c ' (2x - 1)z ' c ' l2v - 1} = 2 *> ri L M1 - (108) 

Performing the external integration first leads to the simple elementary function, so that the final 
integral looks like: 

J « = L du T^W (109) 

It can be easily evaluated by considering its analytical continuation to the whole complex plane 
and choosing the integration contour as displayed in Fig. 0. It reduces to the evaluation of the 
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residues in the n-th order pole in the point u — b = (1 — (2x — l){2y — l))/4y/xxyy and the point 
u = oo: 

f 1 \/u 2 — 1 

J n = 7T Res + Res ^ -f — . (110) 

\u=b u=^j (u _ 6)n v ; 

The first residue provides the contribution proportional to difference of ^-functions, while the 

second one — to the sum. Namely, for the n = 1, 2 which appear in the transformation of the 

one-loop splitting functions we get 

«=<>(u-&)™ 1 1 1 ~ {7a ~ 1 ^ u ~ 1) [0{y -x)- 9(x - y)] , for n = 2 ' 

and 

Reg vgEI ( -1 [g(y ~ *) + ^ - 1/)] , for n = l 

«=°° (m - b) n | - [% - s) + 0(x - j/)] , for n = 2' 

Thus, we finally obtain 

Ji = =^{a^0(y-z)+yx0(a:-y)}, (113) 

y/xxyy 

J 2 = 27r\^-e(y-x) + ^^e(x-y)\. (114) 
[y-x x-y J 

For instance, for the QQ-kernel0 ®®P(z) = (1 + z 2 )/(l — z) we have 

f dz ™P{z) £ ^$-C]{2x - iyC]{2y _1) = I|J_J 2 - & } , (115) 
•>o ^ iVj-d) tt [2yy ]j yy } 

and after little algebra we come to the known non-singlet kernel. 
7.2 and GQ channels. 

We can make use of the same summation formula Eq. ( |107| ) for the mixed channels too. For 
that we have to differentiate both sides of this identity with respect to the first x or the second 
y arguments. Since the anomalous dimensions of the conformal operators differ from the one of 
the usual local operators by the factor 6/j and j/6 for QG and GQ channels, respectively, we 
will introduce these factors by representing the forward evolution kernels as a convolution with 
appropriate functions. The main reason for that is to use the same Eq. ( |107| ), since otherwise 
the evaluation of the infinite series may turn out to be quite complicated. For the QG and GQ 
sectors we define the kernels 

QG P(z) = f 1 dx- f 1 dy QG P(z)5(z-xy), (116) 
JO X Jo 

GQ P(z) = f 1 dxx 2 "- 1 ^ dy GQ P(z)5(z-xy), (117) 
Jo Jo 



3 In this section, for brevity, we omit the colour prefactors in the evolution kernels. 
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and corresponding summation formulae 



E 



w(x\u) 



C](2x - l)z j q+l(2y - 1) = 2 2u ^-^(xxY- 1 * 



'1181 



7T 



x / du(l - u ) 



2\v-l 



XX 



(2x - 1) - (2y - I) J — u 

w yy 



z(l-z 2 ) 



w(x\u + 1) 



E 



(21/ + j)C££(2s - l)|^q(2y - 1) = 2 2 ^ + 1 ^ 



[I - 2a(u)z + z 2 Y+ 2 ' 

(xx) u+ ^ (119) 



7T 6 



x / du(l - u ) 



2\v-l 



yy 



(2y-l)-(2x-l)J—u 



XX 



z{\ - z 2 ) 



[1 -2a{u)z + z 2 Y+ 2 ' 



Since the parity odd channel does not introduce any new features as compared with results avail- 
able in the literature we turn to the vector sector. The integrals appeared in the calculation are 



just the same as before, namely, we get for the ® G P V (z) = (z 2 + (1 — z) 2 ) splitting function 



V , , ^ wix 



dz qg p v (z) x: 

i=i 



M) cn2x ~ 1)z ' cl 



1 ( 

TV 



x -y_ 



\Kvv) 



■h 



i(2y-i) 

I xx Ax — 2y — 1 , xx , . , " 
— z-4 Ji + 4— 2y - l Jo 

2 z/y yy 



(120) 



and for the GQ P 



v , 



1 + (1 - -z) 2 )/^ kernel 



u> a; § 



(3 + i)C/_ 1 (2o;-l)^C?(2y-l) 



'1211 



= - I - ^/?(4y - 2z - 1) Jx + 4xx(2x - 1) Jo) • 
n { Ayy 2 \ yy J 

Discarding the constant term J = | we reproduce the Chase's kernels, however, if the latter is 
taken into account we come to the expressions given in Eqs. fl8ip. 



8 Summary. 

In the present paper we have developed a formalism, which allows to evaluate the n-th order non- 
diagonal part of the exclusive twist-2 evolution kernels starting form the (n— l)-loop approximation 
for the special conformal anomaly matrix, for the flavour singlet channel in Abelian gauge theory. 
The main tools of our analysis are the renormalized anomalous conformal Ward identities and 
commutator constraint obtained from the algebra of the conformal group. Using this approach 
we have derived 0(a 2 ) corrections to singlet anomalous dimensions of the conformal operators. 
This is sufficient to reconstruct the complete ER-BL splitting functions AB V^ l \x,y) in the NLO 
approximation. We have partially checked our calculations by explicit evaluation of the simple 
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two-loop diagrams with fermion vacuum polarization insertions. The formalism we have described 
allows to write in a simple way the solution of the corresponding evolution equations, since the 
corrections to the conformal partial wave expansion are defined by known one-loop evolution 



kernels and ^-functions ||32| . [26| . We have also derived, for the first time, the reduction formulae 
which allow to find the diagonal part of the ER-BL kernels from the corresponding DGLAP 
analogues. This procedure fix unambiguously the form of the diagonal part of the exclusive 
kernels which are exactly diagonal in the physical as well as in non-physical spaces of moments. 

We wish to thank B. Geyer, A.V. Radyushkin, D. Robaschik for valuable correspondence. A.B. 
was supported by the Alexander von Humboldt Foundation and partially by Russian Foundation 
for Fundamental Research, grant N 96-02-17631. 



A Feynman rules. 

In this appendix we summarize Feynman rules which have been used in the calculation of the 
renormalization constants in the main text. Everywhere it is implied that all momenta on the 
lines are incoming. 

For quark and gluon non-local string operators we have 



\ G O v ] 




9af3 1 


[ G o A \ 




ie af3 -+ J 






(QQV 






\q q a 


H 


7+75 J 



f+a;u(k 2 )f+/3;fj,(ki) 



(3 iK 1 k l+ —iK2k2+ -p g— in\k2+~in2ki+ 



(A.l) 
(A.2) 



In the quark sector we have omitted the contributions coming form the phase factor, since they 
are not of relevance for our present consideration. 

The W-operator in parity-even and -odd channels read, respectively, 



\w v ) 


M 


9af3 } 









g — lKiki + ~ _|_ g— «Klfc2+~iK2fcn 



(A.3) 



where fap^itye^ — k a Ep — kpe a is the strength tensor in the momentum space. 

The only modified Feynman rule coming from the differential operator insertion [Ai], required 
to the order we are interested in, is 



V„ = ^(27r) 4 2*<9_<5 (4) (A:i + h + k 3 ). 



(A.4) 
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It should be compared to the usual one for the quark-gluon vertex 

V, = ig ^(27r) 4 6^(k 1 + k 2 + k 3 ). 



(A.5) 



The resulting expression, composed from propagators and vertices we have just derived, should 
be integrated with respect to the momentum of every internal line, i.e. multiplied by the factor 

d d k e 



n 



(2vr) M ' 



(A.6) 



B Generalized ER-BL kernels. 



Recently, there is renewed interest in the deeply virtual Compton scattering [22, 52, 46, 48 . The 



evolution of the corresponding non-forward parton distribution functions |46 

1 



(h'\ F {K U K 2 )\h) 

(h'\ Q O r ( Kl , K2 )\h) 



dx 



dx 



^-IKIX~IK2(C~X) -p ^-iK2X-iKl(£-x) 



G O r (x,(), 
Q O r (x,(), 



(B.l) 
(B.2) 



where the non-local string operators are defined by Eqs. ( |I~3D and the upper (lower) signs corre- 
sponds to the parity-even (odd) channels, is governed by generalized ER-BL evolution equations 

d_ 

dj? 

We perform the calculations and end up with ER-BL type kernels which look like for the parity- 
even sector 



^^- 2 O r (x,C) = -^Jdx'K r (x,x'X)O r (x'X)- 



(B.3) 



QQ K v (x,x',() 
= C F 



(B.4) 



7@ll(^, X — X ) ~\ 7@ll(^ — C) X — X ') + ©111 fa) X - (,X - x') 

ry nn' nn nn' 



® 112 C? x x ) 2 . 1: q(x,2; x x ) 

x [x — Q 



QG K v (x J x',() = 2T F N f 

GQ K v {x, x', C) = C F [{2a/ - C)©Iiifa, x-(,x-x')- x{x - C)©infa, x-(,x-x') 

and for parity odd one 

QQ K A (x, x', C) = QQ K v (x, x', C), (B.5) 
Q G K A (x, x', C) = 2T F N f e\ 12 (x, x-(,x-x'), 

GQ K A {x,x',C) = C F \{2x - ()e° ul (x,x - (,x - x') - x{x - ()e{ u {x,x - (,x - x'j . 



We have used here the generalized G-functions defined by the integral 

r°° da 



6 



lli2-..lr 



\X\ , X 2 , • • • , X n J 



2m 



-a 



II (ax k - 1 + ioy lk . 



(B.6) 



k=l 
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Using the identities 

%k-i x u x ^x^...) = — - — {Q^ jk ..Xx 1 ,x 2 ,x 3 ,. . .) -Glr\ k .( y x 1 ,x 2 ,x 3 ,. . .)} (B.7) 

= — - — {x 2 e™_i jfc ...(xi, x 2 , x 3 , ...)- xi6"_ lfc ...(xi, x 2 , x 3 , . . .)} , 

X\ — x 2 y } 

we can always reduce them to the lowest ones which have the following explicit form 

0?(aO = 0, (B.8) 

0°(x) = S(x), (B.9) 



9(xi)9(-x 2 ) - 9(x 2 )9(-x 1 
B u {x 1 ,x 2 ) = 



(B.10) 

X\ — x 2 

Again the parity even kernels ( p. 4]) possess the off-diagonal pieces in the non-physical sector. 
To cure this problem in appropriate way, so that the forward limit would not be affected, we 
should add the terms 

a -rrv , i >\ , ^ x(C — x) ( (C — 2a;'), for QG-channel _ 

* • = 29 " (I ' X ~ X { AC - .OK - 2x), for GQ-chanuel " ^ 

In the forward limit Q 2 (x) = 5(x), so that these extra terms die out. 

As usual these kernels possess the following DGLAP and ER-BL limits: 

P ( z ) = - K (z, 1, 0), V(x, y) = -K{x, y, 1). (B.12) 
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